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We use a three-body model to investigate breakup reac- 
tions of n Li (n+n+ 9 Li) on a light target. The interaction 
parameters are constrained by known properties of the two- 
body subsystems, the n Li binding energy and fragmentation 
data. The remaining degrees of freedom are discussed. The 
projectile-target interactions are described by phenomeno- 
logical optical potentials. The model predicts dependence 
on beam energy and target, differences between longitudinal 
and transverse momentum distributions and provides abso- 
lute values for all computed differential cross sections. We 
give an almost complete series of observables and compare 
with corresponding measurements. Remarkably good agree- 
ment is obtained. The relative neutron- 9 Li p-wave content is 
about 40%. A p-resonance, consistent with measurements at 
about 0.5 MeV of width about 0.4 MeV, seems to be necessary. 
The widths of the momentum distributions are insensitive to 
target and beam energy with a tendency to increase towards 
lower energies. The transverse momentum distributions are 
broader than the longitudinal due to the diffraction process. 
The absolute values of the cross sections follow the neutron- 
target cross sections and increase strongly for beam energies 
decreasing below 100 MeV/u. 

PACS number(s): 25.60. Gc, 21.60.Gx, 25.60-t, 21.45,+v 



I. INTRODUCTION 

Large efforts are devoted to investigate the properties 
of halo nuclei |0,0] . In particular two- neutron halo nuclei 
have attracted a lot of attention with 6 He (n+n+ 4 He) 
and n Li (n+n+ 9 Li) as the most prominent examples. 
These nuclei are also Borromean three-body systems, 
where all two-body subsystems are unbound The 
two neutrons (the halo) are weakly bound to an ordinary 
nucleus (the core). The halo is spatially extended and 
the two neutrons have a high probability of being out- 
side the core. The core and halo degrees of freedom then 
approximately decouple and three-body models provide 
a good description of such systems. 

The most detailed properties of halo nuclei are ob- 
tained by measurements of fragment momentum distri- 
butions in breakup reactions on stable targets [p|dl7[. 
The projectile energy in these reactions is very large com- 
pared to the energies of the intrinsic motion of the nu- 
cleons in the core, which in turn is much larger than the 
binding energies of the spatially extended halo particles. 
Such high-energy reactions are tempting to describe in 



the sudden approximation where the three-body binding 
is removed without disturbing the motion of the con- 
stituent particles. The three halo particles continue their 
motion independently without any further interaction. 
The resulting momentum distributions then reflects the 
motion in the initial halo bound state of the two neutrons 
and the core. Thus the unchanged initial three-body 
wave function should describe the observed momentum 
distributions. The Fourier transform of the wave func- 
tion indeed approximately reproduce the core momentum 
distribution M|l8-22|. However, the calculated neutron 
distributions are significantly broader than measured. 

Improvements using Glauber theory are possible 
p3| , p4[ , but the neutron distributions within this ap- 
proach are only reported for two- neutron halos in p4| , 
where the two-neutron removal cross sections and the 
related momentum distributions are nicely described. 
However, an alternative is to improve the physically intu- 
itive geometric picture established by the successful sud- 
den approximation. The next step in such a description 
then amounts to instantaneous removal of one halo parti- 
cle (participant) while the remaining two particles (spec- 
tators) continue to interact on the way to detection. This 
modification has only little influence on the core momen- 
tum distribution but affects significantly the distribution 
of the lighter neutron. Several authors suggested that 
this final state interaction between the spectators plays 
an essential role, especially when low- lying resonances 
are present [[l0 19|2f|. Indeed computations then repro- 
duce the measured momentum distributions remarkably 
well when the final state interaction in a consistent cal- 
culation is precisely the same as in the initial three-body 
wave function pl| , p2fl . In addition the model also fairly 
well describes the invariant mass spectra for the two-body 
system consisting of the core and the remaining neutron. 

This successful model requires, however, a better jus- 
tification. The sudden approximation assumes that the 
transition amplitude is proportional to the overlap be- 
tween initial and final state wave functions and the mo- 
mentum distributions are simply proportional to the 
square of this overlap. The participant-target interac- 
tions, implicitly used, are described as the schematic 
black sphere scattering where only absorption is consid- 
ered. The obvious improvement is to use the phenomeno- 
logical optical model to describe the interaction between 
the participant and the target |^,|^|. The qualitative 
improvement is then inclusion of elastic scattering in ad- 
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dition to absorption. 

At the high energies, where these models mostly are 
tested, neutron elastic scatterin g is about three times 
smaller than neutron absorption [B8|. However, the neu- 
tron momentum distribution is much broader for elastic 
scattering than for absorption. The contribution from 
neutron scattering is therefore important. Whether the 
predictions would agree with measurements still remains 
to be studied systematically. The first important step 
is to assume that the optical model describes the inter- 
action between the target and the participant whereas 
the spectators remain undisturbed by the target but still 
continue their motion under the influence of their mutual 
interaction. The complete three-body breakup reaction 
is then described as a sum of these three independent 
contributions. We then neglect processes where two or 
three halo particles interact simultaneously with the tar- 
get 

The three-body model is strictly only valid for struc- 
tureless particles. The effects of extended density distri- 
butions for the constituent particles are small when the 
major part of the three-body wave function is outside 
the radii of all the three particles. However, this require- 
ment is usually not completely fulfilled for halo nuclei. 
The model is then only accurate for the outer part of the 
wave function. Furthermore, the spatial extensions of the 
particles and the target allow geometric configurations 
where more than one halo particle during the collision 
must get close to the target. These configurations should 
be excluded in the process. For one-nucleon halos this 
so-called shadowing effect is known to produce smaller 
absolute cross sections and narrower momentum distri- 
butions |29-31|. For two-neutron halos the shadowing 
effect was included through profile functions in sophisti- 
cated Glauber calculations of three-body fragmentation 
reactions [Q. We account for shadowing by excluding 
the participant wave function inside spheres around the 
two spectators p6j , p7[ . 

The model then consists of an initial three-body halo 
state, reactions caused by the participant-target optical 
potential, a final state with two independent two-body 
systems, i.e. the two spectators and the participant- 
target, and shadowing which excludes the wave function 
within spheres around the spectators. The differential 
cross sections are then products of the participant-target 
cross section and the square of the overlap previously 
used in the pure sudden approximation. The results 
from the successful sudden approximation are essentially 
recovered, but the model now also, via the optical po- 
tential, contains dependence on beam energy and target, 
distinction between longitudinal and transverse momenta 
and also absolute values of all the cross sections. 

The purpose of this paper is to investigate the breakup 
reactions of 11 Li within the three-body model sketched 
above. We shall show systematic computations of a num- 
ber of observables and predict or compare with measure- 
ments. The results are then all correlated as arising from 
the same model with one set of parameters. In section 



II we descibe the model and the method. In section III 
we compare our results with the available experimental 
information and select an interaction with corresponding 
shadowing parameters which best fits the experimental 
data. Section IV presents predictions for a number of 
observables including their energy dependence. Finally 
section V contains a summary and the conclusions. 



II. MODEL AND METHOD 

The spatially extended three-body halo collides with 
a relatively small target at high energy. Then the prob- 
ability that more than one of the constituents interacts 
strongly with the target is small. The differential cross 
section da is then to a good approximation a sum of 
three terms da^ each describing the independent con- 
tribution to the process from the interaction between the 
target and the halo particle i. This is the assumption 
used in the classical formulation for a weakly bound pro- 
jectile We neglect the binding energy of the initial 
three-body bound state compared to the high energy of 
the beam. The reaction is then described as three parti- 
cles independently interacting with the target as if each 
particle were free. 




FIG. 1. Sketch of the reaction and coordinates used. The 
target is labelled by and {i,j,k} label the particles within 
the three-body projectile. Compare with Eq.(jlJ) . 

The process is described as removal of one particle i 
(participant) while the other two particles j and k (spec- 
tators) both survive the reaction undisturbed (see Fig. 
[I]. The participant is either absorbed (stripped) or elasti- 
cally scattered (diffracted) by the target. The final state 
then consists of two independent subsystems, i.e. the 
target plus participant and the two spectators. The in- 
teraction in the final state between the two spectators 
j and k must be the same as in the initial three-body 
bound state. This consistency is previously only reported 
in refs. 

We shall first show how to compute fragmentation 
cross sections of two-neutron halo nuclei on light targets 
where the Coulomb interaction is negligibly small. We 
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then briefly describe how the wave functions in the ini- 
tial and final states are constructed. Finally we give the 
various interactions used in the calculations. 



A. Cross sections 

The masses and the coordinates and their conjugate 
momenta are denoted to, r andp, respectively. The three 
halo particles and the target are labeled by {i,j,k} = 
{1,2,3} and 0, respectively. The relative coordinates 
Tjk,r i, Ti t jk,R and R , see Fig. [l], are defined by 
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The corresponding conjugate momenta are analogously 
denoted by Pj k ,Poi>Pi,jkiPo,ijk and Poi,jfc> Le - 
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These momenta are related in the same way for the final 
states. We use primes to indicate the final state and 
may therefore add primes on all momenta in Eq.(Q). In 
particular p'- k and p' oi are the relative two-body momenta 
and P' = p' 0i j k is the conjugate of R' in the final state. 

We can now easily express the individual momenta in 
terms of three relative momenta and one absolute mo- 
mentum, for example p jk ,Poi>Po,ijk and Pi + Pj + Pk- 
This would be natural for the initial state whereas it 
is more convenient to use the beam momentum and 
p'j k ,p' oi , P' in the final state. We shall work in the center 
of mass system of the projectile where the beam momen- 
tum is — p and p i +p +p k = 0. Using total momentum 
conservation, p' = p Q — p\ — p'^ — p' k , we obtain the in- 
dividual momenta in the final state as 
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The process is described as a participant-target reac- 
tion and two undisturbed spectators, i.e. momentum 
conservation dictates that 



Pj + P'k = Pj +Pk > Po+Pi=Pa+ Pi 



(4) 



In the rest frame of the projectile Pj + p k = p i and 
Pi = Pi jk- Then Eq.(^) can be rewritten as 
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Let us now first consider elastic scattering of the par- 
ticipant i by the light target. We denote the initial 
three-body projectile wave function by \E'(' 7M ) (r^, rij k ), 
where J and M are the total angular momentum and 
its projection on the z-axis. The final state outgoing 
distorted wave functions of the independent participant- 
target and spectator two-body subsystems are respec- 
tively <j& + y, and <Prt + )., y ., , where s ]k = Sj + s k , 

T,j k — Yij + Efe and Si and £j are spin and projection 
quantum numbers of particle i. The corresponding spin 
functions are x Si s', Xs jk T, jk , where we for convenience 
assume a zero spin target. The participant-target inter- 
action is Voi, where we only consider a light target. A 
possible Coulomb interaction can then be neglected. 

The transition amplitude TW of the process, where 
particle i is elastically scattered by the target while par- 
ticles j and k are undisturbed, is most conveniently com- 
puted in the center of mass system of the three-body 
projectile. Using momentum conservation we obtain 



e iP ' R '\V 0i \^ JM) e iPR ) , (6) 



This amplitude basically factorizes into particant-target 
elastic scattering transition amplitude T^? l y, and the 
overlap between initial and final states of the spectators 
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In breakup computations, where only absorption is in- 
cluded in the sudden approximation, the transition am- 
plitude is proportional to the overlap M^ J ^, k s., see 



|21| , |22|] . This previous approximation is therefore still 
valid provided the elastic scattering process can be ne- 
glected. If furthermore the final state interaction between 
the two spectators is neglected the overlap in Eq.(^) re- 
duces to the Fourier transform of the projectile wave 
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function vf»( JM ). This is the approximation used in the 
first attempt to understand these fragmentation reactions 

The differential diffraction (elastic scattering) cross 
section is then given by 



2tt f 6(E' 0i - E, 
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(10) 



where v = po/mo is the velocity of the target seen 
from the projectile rest frame and E oi = an d 
E' Qi = p , Q i /2fj,Qi are the relative energies of particle I 
and the target in the initial and final states. We use 
here the non-relativistic expressions, since the optical 
model is non-relativistic although obtained through a rel- 
ativistic procedure, see the discussion later in the paper. 
The reduced mass of particle % and the target is here 
denoted fiQi. Then energy conservation demands that 
IPoil = |Po*l- The density of final states dirS 1 is given by 



dv 
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(11) 



The differential diffraction cross section in Eq.(10) can 
now also be rewritten in factorized form. When the par- 
ticipant i has spin or 1/2 and the target has spin 0, we 
get, as shown in Appendix A, the expression 
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where the first factor is the differential cross section for 
the participant-target elastic scattering process, 
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the target. We obtain analogously the corresponding dif- 
ferential absorption (stripping) cross section in the same 
factorized form 



d^ ( l(P\p' jk ) 
dP'dp' ]k 

r (0t) 
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where a ab J is the participant-target absorption cross sec- 
tion. The nine-dimensional differential cross section is 
now reduced to six, since the absorped or stripped par- 
ticle inherently is of no interest in the optical model de- 
scription. 

It is conceptually important to realize that the fac- 
torizations in Eqs.(p^) and ( |l5| ) are incomplete, since 
p 0i in the participant-target cross sections is related to 
P' or (p i ju) in Eq. (|l4|) via the momentum conserva- 
tion in Eq.(@). For high energy reactions the factoriza- 
tion is in practice fairly accurate for two reasons. First 
the range of poi-vahies is limited to an interval around 
m ^+m P° determined by the size of the relatively small 
momenta in the motion of the particles within the pro- 
jectile, see Eq.(||). Second the factor arising from the 
participant-target cross section depends only weakly on 
energy for the large beam energies corresponding to this 
rather small range of poi-values. 

Computations of the fragmentation cross sections are 
now essentially reduced to computations of the overlap 
matrix element M^^P E from |2|] (modified to account 
for shadowing by omission of the unwanted cnofiguar- 
tions) and the two-body elastic and absorption cross sec- 
tions determined by the optical model phase shifts. A 
target with zero spin, e.g. an even-even nucleus like 12 C, 
and a neutron with spin 1/2 as the participant particle is 
of particular interest for halo nuclei. We then have P3| 
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and the second factor is the overlap matrix element used 
in the original formulation of the sudden approximation 
for absorption (HQ 

|M s (p iJfc ,p^)| 2 = ^ £ iM^f , (14) 

which is normalized to one, i.e. gives unity after integra- 
tion over Pijk an( l Pjk or equivalently, by use of Eq.(||), 
over P' and p'j k . 

We also consider the other process, where the partic- 
ipant in the sense of the optical model is absorbed by 



where 9 is the angle between p oi and p' oi and 5^ is the 
phase shift of the partial wave I when the total angular 
momentum is j. The functions g and h are given by 
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where Pi is the f th Legendre polynomial. 

The two contributions (stripping and diffraction or ab- 
sorption and scattering) arising from the interaction Voi 
to any measurable cross section is now obtained by in- 
tegration over the unobserved momenta in Eqs.(|l2|) and 
( |l5| ) . The total cross section is given by the sum of both 
contributions and the weight of each of them is directly 
dictated by the optical potential. We shall compute indi- 
vidual as well as relative momentum distributions in the 
final state both along p (longitundinal) and perpendic- 
ular to p (transverse) . 

In addition to momentum distributions we shall also 
compute other observables like the invariant mass spec- 
trum of the two particles in the final state. This mass 
is invariant under Lorentz transformations and there- 
fore independent of coordinate system. In particular, in 
the rest system for the particles j and k the invariant 
mass reduces for the relevant small energies to the non- 
relativistic kinetic energy Ejk = Ej k — p l2 k /2^,jk of this 
two-body system. Then the momentum variable p'j k can 
be substituted by Ejk by use of dEjk — // , <7// , // . 

Another recently investigated observable is the angu- 
lar distribution of the relative momentum between the 
spectators j and k in a coordinate system with the z- 
axis along the center of mass of their total momentum 
in the final state. The decisive variable is then the angle 
between p' jk and p'- + p' k = p l = -p lJk , where the 
latter vector is given in Eq.Q. To compute this angu- 
lar distribution we express the differential cross sections 
as functions of —Pi^k an d p'jk an d integrate over all the 
variables except the angle between these two vectors. 

The necessary integrations require variable changes in 
the description of the final state. The expressions in 
Eq.([|) provide useful identities for this purpose, i.e. 
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(21) 



We have so far only considered the two contributions 
coming from the stripping and the diffraction of particle 
i via the interaction Vqh. However, the reaction in ques- 
tion may also be a result of the other interactions VJy or 
Vbfc. These cross sections are then simply added. Pos- 
sible interference terms are neglected, since they in any 
case are small for spatially extended projectiles. From 
Eqs.([l2|) and ( |l5| ) we expect that the absolute values all 
are of the same order of magnitude as the corresponding 
participant-target cross sections. 



B. Wave functions and shadowing 

The initial three-body wave function ty( JM ) of the pro- 
jectile is obtained by solving the Faddeev equations in 
coordinate space 0]. We use the three sets of hyper- 
spherical coordinates (/?, fii), fli — {on, Cl X i, ^wiii where 
each i is related to a given Jacobi system, see MM . Then 
\fy(JM) j g a gum Q f £ ne three Faddeev components, which 
in turn for each hyperradius p are expanded in a complete 
set of generalized angular functions $>n\p, fli) 

* JM = iE/»wE ®n(p> "«•) > ( 22 ) 

n i=l 

where p~ 5 / 2 is related to the volume element p 5 dpdfli 
with the angular part dVli — sin 2 on cos 2 aidaidQ X idQ y i. 

These angular wave functions satisfy the angular part 
of the three Faddeev equations, i.e. 

h 2 1 

__A 2 $W + Vjk ($V + + $W) 

= ¥-\\n(p)*$ , (23) 
2m p z 

where {i,j, k} is a cyclic permutation of {1, 2, 3}, m is a 
normalization mass, Vjk is the two-body interaction be- 
tween particles j and k and A 2 is the p-independent part 
of the kinetic energy operator. The analytic expressions 
for A 2 and the kinetic energy operator can, for instance, 
be found in Q. 

The radial expansion coefficients f n (p) are obtained 
from the coupled set of "radial" differential equations (J] 

d 2 2m(E-V 3 (p)) | X n (p) | 15 
dp 2 ft 2 p 2 



Xfn(p)= £ \ 2P nn ,— 



fn'(p), (24) 



where Vs is an anticipated three-body potential and the 
functions P and Q are defined as the angular integrals 

3 

Pnn> (p) 

' r 
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1,7=1*' P 



(25) 
(26) 



The continuum wave function 



Ajk+) 

P' lk s jk^jk 



describing 



the two-body spectator system in the final state is ex- 
panded in partial waves [|22f 
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where ( | ) is a Clebsch-Gordon coefficient, Yg mi is the 
spherical harmonic and 3^*(^r) is the angular wave 
function obtained by coupling orbital £ and spin s to 
the total angular momentum and projection j and m. 
The distorted radial wave functions ul Jk . fpl-i., r,-fc) are 
obtained by solving the Schrodinger equation with the 
appropriate two-body potential. When this interaction 
between the two spectators in the final state is neglected 
the expansion in Eq.(p7|) reduces to the usual expansion 
of plane waves in terms of spherical Besscl functions. 

The participant-target interaction is described by the 
optical potential while the spectators remain unaffected. 
The finite extension of the projectile and the target there- 
fore in addition requires exclusion of configurations where 
the spectators pass the target too close to the participant. 
This corresponds to black sphere models describing the 
spectator-target interactions. This so-called shadowing 
strictly requires exclusion of the initial projectile wave 
function in an infinitely long cylinder with the axis along 
the motion of the participant. However, such a cylin- 
der depends on the dynamics of the reaction and omis- 
sion of these events would be technically difficult in large 
scale systematic computations. Instead we approximate 
the shadowing by excluding spheres of the wave function 
where the participant is close to the spectators. This is 
much simpler and has also the appealing feature that ma- 
jor parts of the three-body wave function describing con- 
figurations, where the spatially extended particles are in- 
side the radii of each other, simultaneously are excluded. 
Then the contribution from densities, where halo and 
core nucleons overlap, decreases and consequently the 
possible uncertainty due to the treatment of the Pauli 
principle must be diminished or perhaps completely elim- 
inated || . 

Thus we account for the shadowing effect by substitut- 
ing zero for the initial three-body wave function *( JM ) 
when the distances between participant i and the two 
spectators j and k are smaller than the shadowing pa- 
rameters rnP and r^n) , respectively. Here we indicated 
that one neutron reacts with the target while the core and 
the other neutron are spectators. Instead of this sharp 
cutoff a smooth function, varying from zero at small dis- 
tances to one at large distances, could easily be used to 
eliminate the unwanted geometric configurations. The 
cutoff radii are in any case related to the sizes of target 
and spectators. With the core or nucleon as spectators a 
reasonable parametrization could then be 

r nc = roV 'A 2 / 3 + AT , r nn = ^r^A 2 t /3 + R 2 N , (28) 

where we for simplicity omitted the superscripts ij and 
ik. The mass numbers of target and core are At and A c 
(A c = 9 for 11 Li) and Rn ~ 1 fm is the sharp cutoff ra- 
dius of the nucleon. We choose the parameter tq w 1.26 
fm, which is adjusted to reproduce the few reported ab- 
solute values of two-neutron removal cross sections. This 
somewhat large value of r may reflect a parametrization 



where the range of the nuclear interaction is included. 
The shadowing effect substantially reduces the absolute 
values of the cross sections |2^j2^]. This sensitivity in 
addition to the approximate nature of the treatment of 
shadowing indicates that the predicted absolute values 
are less accurate than desired. On the other hand the 
shapes of the distributions are fairly insensitive to the 
shadowing parameters although the effects are significant 
at the present level of accuracy. 

C. Interaction parameters 

The model described in the previous subsections needs 
specifications of the interactions appropriate for the pro- 
cess under investigation. We shall here focus on the frag- 
mentation of a n Li ( 9 Li + n + n) projectile on a carbon 
target and divide the parameters into five sets, i.e. (i) 
the four two-body interactions between neutron-neutron, 
(ii) neutron- 9 Li, (iii) neutron-carbon and (iv) 9 Li-carbon 
two-body systems and (v) the three-body interaction for 
the 11 Li three-body system. 

The neutron-neutron interaction V nn contains central, 
spin-orbit, tensor and spin-spin interactions, i.e. 

V nn {r) = V c (r) + V ss (r)si ■ s 2 + V T {r)S 12 

~^~Vso(r)£ nn • s nn , (^^) 

where f„„ is the relative orbital angular momentum, Si 2 
is the tensor operator, s\ and s 2 are the spins of the two 
neutrons and s nn = Si + s 2 . We assume gaussian shapes 
for each of the radial potentials V(r). The parameters 
are adjusted to reproduce the experimental low-energy 
neutron-neutron scattering data for s and p-waves. Dif- 
ferent radial shapes and in general more elaborate po- 
tentials have previously been used in the present con- 
text. The results are essentially indistinguishable from 
each other provided the scattering lengths and effective 
ranges remain unchanged. We shall therefore use the 
interactions specified in p4j , where the corresponding s 
and p-wave scattering lengths and effective ranges also 
are given. 

The neutron- 9 Li interaction is not very well known, al- 
though several pieces of information are available. The 
low-energy properties are almost exclusively determined 
by s and p- waves and we shall therefore only include these 
lowest partial waves. The neutron p 3 / 2 -state is bound by 
about 4 MeV in 9 Li and due to Pauli blocking this state 
is unavailable for the valence neutron. We are then left 
with the S1/2 and p^-states. The total neutron angular 
momentum of 1/2 for each of these states then couple 
to the 9 Li-spin of 3/2 resulting in two pairs of spin-split 
states with angular momenta 1 and 2 and parities corre- 
sponding to s and p-states, respectively. 

The interaction in the present three-body calcula- 
tions are determined by the low-energy scattering prop- 
erties, which in turn determine positions and widths of 
resonances and virtual states. These positions are in 
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fact the decisive properties of the interaction and the 
parametrization of the force is in itself rather unimpor- 
tant. We have therefore basically only four parameters of 
physical importance, i.e. the four positions of the s and 
p-states. Furthermore, the fragmentation results are only 
sensitive to the statistically averaged positions of the s 
and p-states, i.e. the spin-splitting, almost unavoidable 
for the strong interaction, does not significantly influ- 
ence the computed differential cross sections, see the next 
section. The relative position of the s and p-states de- 
termines the p 2 -content in the three-body wave function 
of n Li. The s and p-wave content are roughly believed 
to be of the same order as shown in both calculations 
and interpretations of experimental fragmentation data 
0HJl8|j2|^|,^,|||. This amount of p-wave admix- 
ture strongly influences the fragmentation cross sections, 
which in addition also are very sensitive to the binding 
energy or equivalently to the radius of 11 Li. 

Thus we have two crucial parameters (s and p-state po- 
sitions) and two sensitive observables (three-body bind- 
ing and p-wave content). However, this perfect match 
is upset by additional experimental information, espe- 
cially knowledge about a p-resonance in the neutron- 9 Li 
system |T^j3^-^] . As we shall see all three experimental 
constraints can only be reproduced simultaneously by use 
of a three-body interaction, which then is constructed to 
add the missing binding energy in the three-body system. 

We shall therefore first concentrate on the neutron- 
Li two-body system. We assume central, spin-spin and 
spin-orbit terms in the neutron-core interaction, i.e. 



vW(r) = V c W(r)+vW(r)s n -s c 
+V<P(r)t nc -a n , 



(30) 



where £ nc is the relative orbital angular momentum, s„ 
and s c are the intrinsic spins of the neutron and the core. 
As for the neutron-neutron interaction the radial poten- 
tials V(r) tx exp(— r 2 /6 2 ) are gaussians, adjusted inde- 
pendently for each partial wave. 

We choose a large inverse spin-orbit strength to make 
the p 3 /2-interaction sufficiently repulsive to avoid any 
contribution of this partial wave in the three-body wave 
function. Furthermore we shall use a shallow s-wave po- 
tential without bound states. This automatically ex- 
cludes the lowest Pauli forbidden neutron-core s-state 
from the three-body wave function [p|p^]. 

As an important constraint we shall use the "strong 
evidence" for a p-resonance at 538 ± 62 keV with a width 
of 358 ± 23 keV @. We then place the lowest p 1/2 - 
resonance at 0.5 MeV with a width of 0.4 MeV consis- 
tent with the experimental values. We choose rather ar- 
bitrarily the total spin of this state to be 1. The range of 
the gaussian interaction is in this way determined to be 
b = 2.55 fm. A different range produces different widths 
of the resonances. We then continue again somewhat ar- 
bitrarily by placing the other p^-resonance at 0.92 MeV. 
To get a potential at the limit without a three-body po- 
tential we determine the statistically averaged position of 



the two si/2-states to be at 350 keV such that the exper- 
imental n Li binding energy is reproduced. This fixes the 
strength of the s-interaction when we use the same range 
as for the p-interaction. Finally we use the spin-spin in- 
teraction to place the lowest s^-state of total spin 2 at 
an energy of 230 keV. This potential, labeled 77, leads 
to a 30% p-wave content in the 11 Li wave function. We 
could as well have chosen the lowest sm-state with total 
spin 1 at the same energy of 230 keV. The results for 
such a potential would be indistiguishable provided that 
the average energy of the Si/2-states is the same p2| . 

We could also have chosen a different position for the 
highest p^-state. However, a lower value than 0.92 
MeV would place the two px/2 _ states rather close. A 
higher value would require a lower statistically averaged 
position of the Si/2-states, provided the n Li binding en- 
ergy remains unchanged, and the resulting p-wave con- 
tent then becomes unreasonably small (below 30%) in 
disagreement with the fragmentation data [fio |l4 , 24[| . 

We can now increase the p-wave content by increas- 
ing the statistically averaged position of the s^-states 
while the parameters for the p-waves remain unchanged. 
The resulting undcrbinding of 11 Li must then be com- 
pensated by an attractive three-body force. In this way 
we construct potential I with 40% p-wave content. Again 
we could have chosen a different position for the highest 
Pi/2-state while maintaining a reasonable p-wave content 
between 35% and 40%. However, a lower value would 

TABLE I. Parameters for various neutron- 9 Li interactions. 
The form is given in Eq.(]3C|) and the radial shapes are all gaus- 
sians, exp(— r 2 /b 2 ), with ranges b — 2.55 fm and strengths de- 
noted by Vc , VjP and Vii' with £ — 0,1 for s and p-waves, 
respectively. The strength parameter for the corresponding 
three-body interaction, V3 exp( — p' 2 /b 2 ) with 63 = 2.50 fm, 
are also given. The lower part of the table contains p-wave 
content in % in the n Li wave function and the four energies 
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then require a lower statistically averaged position of the 
Si/2-states and a repulsive three-body force in order to 
get the correct 11 Li binding energy. A higher value would 
require a rather high statistically averaged position of the 
sj/2-states, i.e. close or even above the statistical average 
of the pi/2-states. The freedom in choosing the positions 
of the resonances and virtual states is in fact limited. 

The parameters of the potentials I and II are given in 
Table |I along with the p-wave content and the energies 
of the resonances and virtual states. The statistically 
averaged position of the two virtual s-states in 10 Li is 
higher for potential I than for potential II resulting in 
the correspondingly larger p-wave content. 

In Table | we also for comparison give the parameters 
corresponding to the neutron- 9 Li interaction (called III) 
used in fh]| , f4l| ] . The main difference compared to poten- 
tials / and II is that the lowest p-resonance is placed 
at 0.75 MeV with a width of 0.87 MeV. Maintaining the 
p-interaction and increasing the p-wave content to about 
40% would require a statistically averaged position of the 
s-states above the corresponding average position of the 
p-states. This is the reason for the choice of the slightly 
lower value of 0.5 MeV for the lowest p-state. Potential I 
is a rather good starting point. To investigate the effects 
of spin-splitting we also constructed the potentials IV 
and V, where the statistically averaged positions of the 
s and p-waves are maintained. Then the p-wave content 
and the 11 Li binding energy is also unchanged. 

In conclusion, the neutron- 9 Li potential is already 
rather severely constraint by existing data. The lowest 
Pi/2-state is experimentally determined to be around 0.5 
MeV. The ground state is an s-statc. The distance be- 
tween the statistically averaged positions for the s 1 / 2 an d 
Pi/2-states is determined by the requirement of about 
40% p-wave content of the neutron-core relative motion 
in the total 11 Li wave function. For example maintaining 
the 11 Li binding energy at around 300 keV, the p-wave 
content around 40% and the three-body force in potential 
L a P1/2 average energy at around 0.65 MeV would re- 
quire the Sim average energy at roughly the same value. 
If the lowest p^-state is at 0.5 MeV the highest pi/2~ 
state has to be above 0.72 MeV, because the pi/2 aver- 
age energy otherwise would fall below the S1/2 average 
energy. In the same way, a low s^-state close to zero 
requires the highest s^-state below 1.7 MeV to avoid 
the same inversion. 

The three-body interaction is chosen as V3 exp(— p 2 /6|), 
where 63 = 2.50 fm and the strength V3 is adjusted to 
reproduce the measured two-neutron separation energy 
for 11 Li, i.e. 295 ± 35 keV. This additional force is nec- 
essary, since the two-body interactions reproducing all 
low-energy scattering phase shifts lead to an undcrbind- 
ing of the three-body system. For 6 He this deficiency 
amounts to around 0.5 MeV fUJH]. The idea is that 
the three-body force should account for the polarization 
of the particles beyond that described by the two-body 
interactions. Then all three particles must be close to 



produce this additional polarization or modification of 
the intrinsic structure of the composite particles. Thus 
the three-body force must be of short range. 

The importance of the three-body interaction is per- 
haps seen most clearly in three-body continuum calcu- 
lations where the resonance structure of the two-body 
subsystems probably is decisive. The three-body force is 
designed to give the correct three-body binding energy 
while the two-body interactions remain unchanged still 
reproducing the two-body structure. With the correct 
two-body interactions the computed three-body contin- 
uum structure is much more reliable fl43| . 

For the neutron- carbon interaction we use non- 
relativistic optical potentials obtained from relativistic 
potentials through a reduction of the Dirac equation into 
a Schrodinger-like equation j44j|. These phenomenolog- 
ical potentials in the Schrodinger equation produce the 
same scattering data as obtained by_use of the relativistic 
potentials in the Dirac equation 
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FIG. 2. The total (solid), elastic (dashed) and absorption 
(short-dashed) cross sections a for neutron- 12 C reactions as 
functions of the neutron laboratory kinetic energy Ti a t- The 
inset focuses on the region of the minimum. The open and 
filled circles are experimental points from [[|6],0] , respectively. 
The curves are obtained in an optical model computation as 
decribed in using the parameters in J2j]]. We include 35 
partial waves in the calculations. 

In particular for a carbon target the interaction used is 
the parameterization EDAI-C12 |28| valid for a range of 
neutron kinetic energies from 29 to 1040 MeV. The result- 
ing neutron- 12 C cross sections are shown in Fig. |2[ For 
energies above 150 MeV the cross sections are relatively 
weak functions of energy and absorption is much more 
likely than diffraction. The factorization of the fragmen- 
tation cross sections in Eqs.([l2]) and ( fi"5| ) is then rather 
accurate at these high energies. 

At energies below 150 MeV the neutron-target cross 
sections increase dramatically. The factorization is then 
less accurate, but still a useful approximation. Further- 
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more, diffraction quickly dominates over absorption. The 
absolute cross sections for fragmentation of n Li can then 
be expected to increase with decreasing beam energy. 
The shapes of the distributions also should change from 
narrow (absorption) to broad (elastic scattering) with de- 
creasing energy. At intermediate energies we can expect 
a mixture exhibiting a narrow distribution with a rel- 
atively large broad background. These predictions are 
simple consequences of the model, where the fragmen- 
tation cross sections essentially are proportional to the 
neutron-carbon cross sections. 

The carbon- 9 Li interaction is needed to estimate the 
cross sections both when the 9 Li-core is destroyed (ab- 
sorped) and the two neutrons detected and when the 
core is scattered (diffracted) on the target. These cross 
sections are obtained by a two-body computation of the 
carbon- 9 Li reactions using the simple almost schematic 
optical potential defined in |Q and described in [Q . In 
this model the beam energy dependence is introduced 
through the proton-proton and neutron-proton cross sec- 
tions. Experimental data of these nucleon-nucleon cross 
sections can be found in pCj ], 

In the calculation of the 9 Li- 12 C cross section we as- 
sume spin zero for 9 Li. We need around 150 partial waves 
to get convergence. At a beam energy of 280 MeV/u the 
computed elastic cross section is 419 mb and the com- 
puted absorption cross section is 795 mb. This gives a 
total of 1214 mb consistent with the results in J5l|]. These 
values are dramatically reduced by the shadowing effect. 

III. COMPARISON WITH EXPERIMENTS 

The interaction parameters discussed in the previous 
section determine the structure of the three-body pro- 
jectile. The positions of the resonances and the p-wave 
content are the crucial quantities. The reactions are de- 
scribed by the participant-target interactions and the 
cutoff radii taking the shadowing effect into account. 
We previously used the cutoff radii r nc = r nn = 3 fm 
JE7J , which is consistent with the fragmentation data for 
°He on carbon p6fl . For n Li fragmentation the neutron- 
neutron shadowing parameter should remain the same, 
since the physical origin is unchanged. Therefore we 
use r nn = 3 fm while we still maintain the neutron- 
9 Li shadowing radius r nc as a parameter for adjustments 
within rather narrow limits around the value obtained 
from Eq. (|28|) . The radius of 9 Li is 0.9 fm larger than the 
radius of ^le and the value of r nc is then expected to 
be about 4 fm for 11 Li. In this section we shall compute 
different types of observables as discussed in the previous 
section and use the experimental data to select a promis- 
ing set of interactions and shadowing parameters. We 
shall only display results for potentials I, II and III and 
omit the curves for potentials IV and V, which in all cases 
barely can be distinguished from those of I. 

The two-neutron removal cross section <r_2n is known 



experimentally for 11 Li fragmentation on a carbon target 
at 280 MeV/u All events, where 9 Li is detected after 
the fragmentation, contribute to this cross section. Ac- 
cording to fl4|| three different reaction mechanisms lead 
to such a halo breakup. First electromagnetic excitation 
of the halo state into the continuum followed by decay of 
these excited states. We shall neglect this process, since 
it is expected to give a rather small contribution (less 
than 10 mb) for a light target like carbon. 

The second mechanism is the stripping of one of the 
halo neutrons by the target. This process is also called 
absorption and the corresponding cross section is denoted 
by cr_ 2n . It is fully described in our model via the imagi- 
nary part of the neutron-target optical potential. On the 
other hand we neglect processes where both neutrons are 
absorbed by the target and no neutrons appear in the 
final state. This contribution to the cross section is ex- 
pected to be small, since the large spatial extension of 
the projectile diminishes these simultaneous reactions. 

TABLE II. Core diffraction (er^^), two-neutron removal 
cross sections (cr_2n) equal to the sum of neutron diffraction 
(°"-2n) an d neutron stripping (o"^2n) computed for fragmen- 
tation of lx Li at an energy of 280 MeV/u. The lowest part of 
the table is for an aluminum target whereas everything else 
refers to a carbon target. The cross sections are in millibarns. 
The potentials and shadowing parameters are from Table [l| 
The experimental data are from jl4| and |L2| for carbon and 
aluminum, respectively. Potentials IV and V produce the 
same values as potential I. 
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The third mechanism is the diffraction of one of the 
halo neutrons. We also refer to this process as elastic 
scattering and we denote the corresponding cross section 
by cr? 2n - For weakly bound systems and not too small 
beam energy the main contribution to cr? 2n comes from 
processes where the neutron after scattering on the target 
ceases to interact with the remaining two halo particles. 
This is precisely our model, where the final state is de- 
scribed as two independent two-body systems and the 
interaction between them is ignored, see Fig. |l|. 

In addition to these three mechanisms there must be 
a contribution (crSj where the core is scattered by the 
target. In Table O we give these computed cross sec- 
tions for the potentials in Table | and different sets of 
shadowing parameters. It is remarkable that the cross 
sections calculated without shadowing are virtually inde- 
pendent of the potentials, and furthermore clearly much 
larger than observed. The computed values are reduced 
by roughly a factor of two by using the shadowing pa- 
rameters r nc — 3.5 fm and r nn = 3 fm. Then by neglect- 
ing the two-neutron absorption (not computed) and the 
core diffraction processes the experimental numbers are 
reproduced within the error bars for the total as well as 
for the stripping and diffraction cross sections. Therefore 
we should underestimate the experimental values and use 
the larger shadowing parameters r nc = 4 fm and r nn = 3 
consistent with the parametrization in Eq. (p8|) . 
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FIG. 3. Radial neutron momentum distribution after frag- 
mentation of n Li on a carbon target at 280 MeV/u. Core 
diffraction is not included. The experimental data from 
is compared with calculations using the potentials I (solid), 
II (short-dashed) and III (long-dashed) in Table |. The shad- 
owing parameters are given on the figure. The optical model 
parameters are from Eg . 

However, the computed core diffraction contributions 
in Table |n] are not negligible, but for comparison with ex- 
perimental values they should be reduced, since (i) a frac- 
tion corresponds to elastically scattered 11 Li-particles, 
(ii) the large transverse momenta of 9 Li precisely ob- 
tained through the scattering process to a large extent are 
experimentally excluded due to limitations in the large 
angle acceptance, (iii) the shadowing by the core is prob- 



ably larger than the shadowing by a neutron correspond- 
ing perhaps rather to r nc — 5 fm and r nn = 4. Further- 
more, the calculations rely on rather uncertain optical 
model parameters. The effect is almost entirely confined 
to the absolute values while the shapes of the distribu- 
tions only are marginally influenced. Using the complete 
factorization approximation we estimate the size of this 
cross section to be about 60% of the values of cr_2 n m 
Table [n| If necessary it can be added at the appropriate 
places, but we shall in the remaining part of this paper 
not include this process unless explicitly mentioned. 

The radial neutron momentum distribution after neu- 
tron removal in fragmentation of a 280 MeV/u Li- 
projectile on carbon is measured |1C| ], The variable is p r , 
where p 2 r = p%+Py is expressed in terms of the projections 
p x and p v of the neutron momentum along the two direc- 
tions perpendicular to the beam direction chosen as the 
z-axis. In Fig. [}] we compare our calculations (suitably 
scaled) for different potentials with the measured mo- 
mentum distributions in arbitrary units. Without shad- 
owing, shown in part (a), potentials / and III give too 
broad and too narrow distributions, respectively, whereas 
potential II with 30% p-wave content reproduces the 
data. This result is consistent with previous computa- 
tions without shadowing and diffraction |2l],|22|]. With 
shadowing, shown in parts (b) and (c), we obtain as ex- 
pected narrower distributions. Variation from r nc = 3.5 
fm to 4 fm is hardly visible showing a very weak depen- 
dence on the r nc shadowing parameter. Now potential 
III gives a too narrow neutron momentum distributions 
while potentials I and II both reproduce the data. 

The invariant mass spectrum of 1 Li is independent of 
the coordinate system and in the rest system of 10 Li the 
invariant mass is (after subtraction of the rest masses) 
equal to the total kinetic energy E nc of the neutron- 9 Li 
system Jl4|,|34| . In Fig. ||| we compare calculations with 
experimental data obtained after fragmentation of 11 Li 
p4| . The experimental distribution is given in absolute 
numbers and we compare directly without any arbitrary 
scaling. We only include processes where the spectators 
form the detected 10 Li, i.e. we neglect those neutron-core 
combinations where the neutron or the core have been 
through the scattering process. The latter contributions 
are relatively small, i.e. the core diffraction and half of 
the neutron diffraction cross sections. Furthermore the 
corresponding contributions from these participants are 
probably not fully included in the measurement. 

The calculations without shadowing in panel (a) pro- 
duce for all potentials distributions with maxima much 
higher than observed. With shadowing the maxima are 
reduced as shown in panels (b) and (c), which again are 
rather similar, but still potentials 77 and 1/7 both give 
too high peaks. Clearly the best comparison is obtained 
for potential I and especially with the shadowing in panel 
(c). The discrepancies appear as a slightly too low-lying 
peak energy and a slightly too narrow peak. However, 
the response function of the detector system is not ac- 
counted for and a more accurate comparison would shift 
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the peak towards higher values, decrease the maximum 
value and broaden the peak Jl4| . Thus all defiencies may 
be improved in this way. In any case the behavior of 
the invariant mass spectrum is strongly influenced by 
the properties of the low- lying resonances. For exam- 
ple a narrow p-resonance would produce a shoulder or a 
peak in the distribution at the position of the resonance. 
A broader resonance would only show up as a larger cross 
section at the corresponding position. Thus a more de- 
tailed and more accurate comparison between computed 
and measured distributions would provide information 
about the neutron-core resonance structure. 
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FIG. 4. Invariant mass spectrum of 10 Li after fragmenta- 
tion of lx Li on a carbon target at 280 MeV/u. Only neutron 
and core spectators forming 10 Li are included. The experi- 
mental data are from [jl4j. The calculated curves in (a), (b 
and (c) are obtained with the same parameters as in Fig. 
The maxima for the computed curves in panel (a) (outside 
the figure) are around 500 mb/MeV, 650 mb/MeV and 1000 
mb/MeV for potentials I, II and III, respectively. In panel 
(b) and (c) the maxima for potential II are 430 mb/MeV 
and 400 mb/MeV, respectively. For potential III we obtain 
instead 650 mb/MeV and 600 mb/MeV. 

The radial neutron momentum distribution after core 
breakup of 11 Li on a carbon target at 280 MeV/u is 
also known experimentally for collisions where 9 Li is de- 
stroyed |rlj| . The neutrons are detected in coincidence 
with a charged fragment different from 9 Li. The contri- 
bution of neutrons coming from the core is eliminated in 
the experiment by subtraction of the corresponding neu- 
tron momentum distribution obtained from the measured 
9 Li-carbon reaction. 

This distribution arises only from absorption, i.e. we 
must compute the cross section in Eq.(|lE|), where the 
9 Li-core is destroyed by the target while the neutrons 
continue unaffected. For simplicity we assume complete 
factorization such that the 9 Li-carbon interaction only is 
necessary to provide the absolute scale through the ab- 
sorption cross section. This is a very good approximation 
for large beam energies. The momentum distribution 
is then given by Eq.(|l4|) and therefore computed as de- 
scribed in |p2| . The absolute cross section is obtained af- 
terwards by multiplication of the absorption cross section 



computed with the optical model parameters in |t8 49 
to be 795 mb for a beam energy of 280 MeV/u. 

Both from Table || and Figs. || and || we found 
good agreement between theory and experiment with the 
choice r nc — 4 fm and r nn — 3 fm for the shadowing pa- 
rameters. For core breakup reactions, where only the 
neutron-core shadowing parameter is relevant, we there- 
fore use r nc = 4 fm in the computation of the neutron 
momentum distribution shown in Fig. |^. The agreement 
between theory and experiment is remarkably good in 
view of the simplicity of the model, where the interactions 
between the halo neutrons and all the fragments from the 
core destruction have been neglected. In contrast to the 
other observables this reaction produces fragments with 
approximately the same velocity as the halo neutrons and 
final state interactions beside that of the two neutrons 
could be significant. The two spectator neutrons could 
also be disturbed during such violent reactions. 
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FIG. 5. Radial neutron momentum distribution after core 
breakup fragmentation of n Li on a carbon target at 280 
MeV/u. The calculations are for the interactions in Table 
I and a neutron-core shadowing parameter of 4 fm. The 
Li-carbon interactio n pro viding the absolute values is the 
optical model from ^8|4s|. The experimental data Jllj in 
arbitrary units are scaled to match the calculations. 

The larger p-wave content arising for potential / gives 
a broader distribution closer to the data than the poten- 
tials II and 177. The comparison in Fig. || is then also 
supporting the choice of potential I and a neutron-core 
shadowing parameter of 4 fm in the description of the 
fragmentation reactions of 11 Li on a light target. 

The transverse 9 Li- core momentum distribution for 
breakup of n Li on a carbon target at 280 MeV/u is com- 
puted and shown in Fig. Jp| for the potentials and shadow- 
ing parameters in Table BT Detailed comparison with the 
data |l2|,|l5| requires folding of the computed curves with 
the experimental beam profile resulting in a few MeV 
broader curves p0[ . Due to the large experimental er- 
rors this distribution is not very helpful in constraining 
the potentials and shadowing parameters. However, we 



11 



can still conclude that the agreement between theory and 
experiment is satisfactory for potential /. 

In Fig. ^ we also show the experimental longitudinal 
core momentum distribution on an aluminum target for 
two different beam energies. In the next section we shall 
show that the difference between longitudinal and trans- 
verse momentum distributions is rather small and the 
widths of the distributions are almost independent of the 
beam energy. The difference between the transverse and 
longitudinal data in Fig. || can be explained by the use 
of a different target in the experiments. The aluminum 
radius is almost one fermi larger than the carbon radius, 
the optical model parameters are different and more im- 
portant the shadowing parameters must be larger for alu- 
minum with the resulting narrower momentum distribu- 
tions as seen in Fig. |[ Still the tail of the distribu- 
tion is not reproduced. The theoretical prediction of the 
two-neutron removal cross sections for fragmentation of 
11 Li on an aluminum target are shown in Table 0. The 
parameters (r nc ,r nn ) — (5 fm, 4 fm) are consistent with 
Eq. ( p8| ) and the corresponding two-neutron removal cross 
section is then expected to be c_2n ~ 350 mb. 
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FIG. 6. Momentum distributions for the 9 Li-core after 
fragmentation of 11 Li. Core diffraction is not included. The 
left hand side is the transverse distribution for a beam energy 
of 280 MeV/u on a carbon target and the experimental data 
is from jUI (filled squares) and @ (filled triangles). The cal- 
culations are for the potentials in Table | with (r nc , r n „) = (4 
fm, 3 fm). The right hand side shows the experimental longi- 
tudinal distribution for a beam energy of 468 MeV/u (filled 
circles) and 648 MeV/u (open circles) on an aluminum target 
[ pTs| . The calculations are for an energy of 280 MeV/u for po- 
tential / and the three sets of shadowing parameters specified 
on the figure. The optical model parameters are from [ p8{ . 
All the computed distributions are scaled to the experiment. 

The conclusion of this section is that potential I with 
the shadowing parameters (r nc , r nn ) — (4 fm, 3 fm) for a 
carbon target gives an excellent agreement between the- 
ory and experiment for the observables discussed. This 
is a strong justification for the model and the method. 



IV. BEAM ENERGY DEPENDENCE 

We have now established a model with a set of param- 
eters successfully reproducing a variety of experimental 



data, i.e. potential I, r, n 



4 fm, 



= 3 fm and op- 
tical model parameters from |2^|49|]. In the remaining 
part of this paper we shall only use these parameters and 
explore the consequences of the model for a number of 
observables. Particularly we shall in the following con- 
centrate on the predicted energy dependence of various 
quantities in fragmentation reactions of n Li on carbon. 
The model only has a dependence on the beam energy 
through the interaction between the participant and the 
target. This interaction is described by phenomenological 
optical models, which give absorption and elastic scatter- 
ing cross sections as functions of particle energy, see Fig . 
||. These cross sections are decisive factors in Eqs.(|l~2|) 
and ([l5]) and two-neutron removal cross sections must 
show the same energy dependence. 
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FIG. 7. Two-neutron removal cross sections for 11 Li frag- 
mentation on carbon (thick curves) and aluminum (thin 
curves) as functions of the beam energy. The cross sec- 
tions from neutron stripping a^. 2n (short-dashed) and neu- 
tron diffraction a^2n (long-dashed) are shown along with 
the SLllll (7 — 2n (solid). The shadowing parameters are 
(r nc ,r nn ) — (4,3) fm for carbon and (5,4) fm for aluminum, 
the interactions are potential I from Table @ and the optical 
model is specified in |28^9| . The experimental data are from 
flftipt for increasing energies for carbon and from |l^] for 
aluminum. For |^,^| only the total cross section is given. 

This is indeed seen in Fig. ^ where we show two- 
neutron removal cross sections as functions of the beam 
energy. These calculations do not include two-neutron 
absorption processes and processes where the core inter- 
acts with the target. The cross sections are clearly gov- 
erned by the behavior of the neutron-carbon cross sec- 
tion, with a minimum at a beam energy of around 250 
MeV/u. For larger energies we observe smooth increases 
towards a flat region. For smaller energies the cross sec- 
tions increase rather dramatically. The prediction is an 
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increase by about 70% when carbon and aluminum are 
interchanged as target. The computed curves underesti- 
mate the latest experimental points [Q , as expected due 
to the neglect of core diffraction. On the other hand the 
older data are far below the calculations. However, at 
30 MeV/u total two-neutron removal cross sections are 
measured for targets of beryllium 0.47 ±0.10 b and nickel 
1.3 ± 0.4 b 0. This is in agreement with the computed 
increase towards smaller beam energies. 

The two-neutron removal cross sections can be sepa- 
rated into a number of differential cross sections. Wc 
shall discuss the momentum distributions of the halo par- 
ticles both relative to each other and individually with 
respect to the center of mass of the projectile. In addi- 
tion we shall discuss the invariant mass of 11 Li and the 
angular correlation of the emitted neutrons. These ob- 
servables are in most cases not experimentally available 
and our results are therefore model predictions. 
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FIG. 8. Longitudinal (thick) and transverse (thin) neutron 
momentum distributions in coincidence with 9 Li after frag- 
mentation of n Li on carbon at 280 MeV/u computed in the 
rest frame of the projectile. The beam momentum is the ref- 
erence direction. The short-dashed and long-dashed curves 
are the contributions to the total (solid) from diffraction and 
stripping of the participant neutron. Core diffraction is shown 
separately as the filled circles. The inset shows the total trans- 
verse neutron momentum distribution for the same reaction 
for different beam energies in the same unit. The standard 
set of parameters in Fig. |^ is used. 

In Fig. ^| we show the neutron momentum distribu- 
tions after fragmentation of n Li on carbon. When the 
participant neutron is absorbed in the stripping process, 
the detected neutron must be a spectator and then there 
is no difference between longitudinal and transverse mo- 
mentum distributions in the rest frame of the projectile. 
When the participant neutron is scattered by the target 
(diffraction) it receives additional momentum perpendic- 
ular to the direction of the beam. This process therefore 
contributes with a broader momentum distribution in the 



transverse than in the longitudinal direction. The tail in 
the total neutron momentum distribution is then more 
pronounced for the transverse than for the longitudinal 
distribution. On the other hand different tails do not 
necessarily imply that the full width at half maximum 
(FWHM) of the distributions also differs. The relative 
size of the stripping and diffraction contributions reflects 
the size of the corresponding cross sections in Fig. [?]. The 
contribution to the neutron momentum distribution from 
the core diffraction (dotted line in the figure) is the same 
in the longitudinal and transverse directions, as expected 
because both neutrons are spectators. This contribution 
of about 30 mb is not added in the figure to the total, 
but the total width would only be marginally influenced. 

The inset in Fig. || shows the transverse neutron mo- 
mentum distribution for different beam energies. The 
shape of the distributions is essentially independent 
of the energy due to the approximate factorization in 
Eqs.(^2|) and ( |l5[ ) at these fairly high energies, see Table 
III . The computed FWHM decreases slightly with an en- 
ergy increase from 50 MeV/u to 280 MeV/u and remains 
then essentially constant at higher energies. In contrast, 
the maximum or peak value for the momentum distribu- 
tions changes considerably with the beam energy. The 
behavior of the peak values, i.e. a sharp decrease, the 
passing of a minimum around 250 MeV/u followed by a 
smooth increase is as expected similar to the variation 
shown in Fig. 0. 



TABLE III. The maximum values in mb/MeV/c 
(one-dimensional) and mb/MeV 2 /c 2 (radial) and the FWHM 
in MeV/c as functions of beam energy in MeV/u for various 
momentum distributions computed for fragmentation of 11 Li 
on carbon in the rest frame of the projectile. Core diffraction 
is not included. The indices _L, || and r indicate transverse, 
longitudinal and radial distributions, respectively. For the ra- 
dial distributions the FWHM is the width at half maximum 
multiplied by two. For the invariant mass spectrum we give 
the maximum value in mb/MeV and the width in MeV. 
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Instead of referring the neutron momentum distribu- 
tion to the rest system of the projectile as in Fig. [§[ we 
could refer it to the rest system of the spectator neutron- 
9 Li system. Actually in this frame the neutron momen- 
tum is the relative neutron- 9 Li momentum. This momen- 
tum distribution is shown in Fig. ^, and it is identical to 
those of Fig. || for an infinitely heavy core. In Fig. |9| we 
only include the dominating term arising from the specta- 
tor neutron and the longitudinal and transverse distribu- 
tions are therefore identical. The neglected contribution 
is about half of the diffraction part. The FWHM and the 
peak values for the curves in the inset are given in Table 
Bj|. The FWHM is 37 MeV/c for all the energies coincid- 
ing with the energy independent width of the stripping 
part in Fig. ||. 

For comparison the neutron-neutron relative momen- 
tum distribution computed using the complete factoriza- 
tion approximation is also shown in Fig. |^. The shape is 
in this approximation energy independent and the scale is 
determined by the carbon- 9 Li absorption (core destruc- 
tion) cross section. For diffraction (core survival) we 
should multiply it by about a factor 0.53 for a beam 
energy of 280 MeV/u. 
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FIG. 9. The neutron- 9 Li relative momentum distribution 
for the cases in Fig. |§| The relatively small contributions 
from the neutron and core participants are not included and 
the longitudinal and transverse momentum distributions are 
therefore identical. The neutron-neutron relative momentum 
distribution computed for core destruction using the complete 
factorization approximation is also shown for an energy of 
280 MeV/u. The carbon- 9 Li absorption and diffraction cross 
sections are 795 mb and 419 mb, respectively. 

The two-dimensional radial momentum distributions 
are often used to increase the number of observed events. 
The variable is then p r (p^ = p\ + Py) and integration 
of this momentum distribution over p x (or p v ) gives the 
transverse momentum. The results are shown in Fig. [l^ 
for the cases in Fig. ^. The widths and the peak values 
for the curves in the inset are given in Table III . We find 



the same qualitative behavior as for the one-dimensional 
distributions. This variation as well as the relative size 
of the stripping and diffraction contributions are again 
consistent with the result in Fig. @. 
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FIG. 10. The 
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neutron- 9 Li and the neutron-neutron relative momenta for 
the cases in Fig. []. 

The radial neutron momentum distribution analogous 
to Fig. ^ is almost indistinguishable from the results in 
Fig. nq. The cross section for the neutron-neutron mo- 
mentum distribution corresponds to core destruction for 
a beam energy of 280 MeV/u. The core survival process 
is obtained by multiplication with the factor 0.53. The 
widths are energy independent in this approximation. 
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FIG. 11. The 9 Li momentum distribution for the cases in 
Fig. ^. The transverse and longitudinal contributions from 
the diffracted core are shown separately, but not included 
in the total cross section. The longitudinal and transverse 
momentum distributions are therefore identical. The inset 
shows the distributions in the same unit for different beam 
energies. 
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The neutron momentum distribution is narrower than 
that of the 9 Li core due to the final state interaction. We 
show in Fig. |ll| the computed 9 Li momentum distribu- 
tions for the same cases as in Fig. |[ The transverse 
distribution from the diffracted core is small and very 
broad due to the diffraction process. The longitudinal 
distribution is as usual narrower. These contributions 
add about 30 mb to the total cross section while chang- 
ing only marginally the shape of the total distribution. 
They are not added in the figure where the total distribu- 
tion then only includes contributions from the participant 
neutrons. The displayed longitudinal and transverse mo- 
mentum distributions are therefore identical, since the 
difference between them is due to the diffraction process. 
The inset shows the core momentum distributions for 
different beam energies. The computed widths and the 
peak values are given in Table |ff]|. The behavior is again 
a reflection of the results in Fig. |2|. 

The momentum distributions of 9 Li are not far from 
those found in the simplest approximation described by 
the Fourier transform of the initial three-body wave func- 
tion. However, the neutron momentum distributions are 
strongly influenced by the final state interaction. Instead 
the momentum distribution of the center of mass of 10 Li 
shown in Fig. |l2| reveals direct information about the 
neutron momentum distribution in the initial three-body 
system. The process is removal (stripping or diffraction) 
of one halo neutron by the target with the remaining 10 Li 
system as the two spectators. In the center of mass of 
the three-body projectile the momentum distributions of 
10 Li and the participant neutron are then identical. 
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FIG. 12. The momentum distribution for the Li center of 
mass for the cases in Fig. ^. The relatively small contributions 
from the neutron and core participants are not included and 
the longitudinal and transverse momentum distributions are 
therefore identical. The inset shows the distributions in the 
same unit for different beam energies. 

This observable is insensitive to the final state inter- 
action where the opposite erroneously was postulated in 



p2| . On the other hand the sensitivity to the shadow- 
ing parameters is large Jl6|,^6| . Experimental data could 
then be very useful to check the validity of the shad- 
owing parameters extracted in the previous section from 
the comparison between computed momentum distribu- 
tions and available experimental data. The final state 
momentum of the participant neutron does not enter in 
the measured momentum and the longitudinal and trans- 
verse 10 Li momentum distribution are therefore identi- 
cal. This can also be understood from the fact that the 
initial three-body momentum of the participant neutron 
does not have a preferred direction. The inset of Fig. 
|l2| again shows the variation of the distribution with the 
beam energy. The behavior is the same as discussed in 
connection with the previous figures. The related key 



numbers are given in Table III 
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FIG. 13. The distribution of the relative angle 6 in the rest 
frame of the initial three-body system. The angle 8 is defined 
as the angle between the center of mass momentum of 10 Li 
and the relative neutron- 9 Li momentum. The reactions are 
the same as in Fig. ^. The relatively small contributions from 
the neutron and core participants are not included. The inset 
shows the distributions for different beam energies. 

The momentum distributions were recently supple- 
mented by an angular correlation measurement, where 
the observable is the angular distribution of the relative 
momentum between the detected neutron and the core in 
a coordinate system with the z-axis along the center of 
mass momentum of the neutron-plus-core spectator sys- 
tem jL7|. This observable is shown in Fig. [l3| for our 
standard case of n Li fragmentation on carbon. In the 
computation we assumed that 9 is constructed as the 
angle between the momentum of the center of mass of 
the 10 Li spectator system and the relative momentum 
between 9 Li and the spectator neutron. The small con- 
tributions, where the neutron or the core in 10 Li are the 
scattered participants, are not included in this figure, be- 
cause they are small and furthermore almost completely 
excluded in the experiment. An estimate of the shape 
and size of one of the neglected contributions can be 
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found in p7j , where both 6 He and n Li fragmentation 
are discussed. 



The computations as always involve Eqs.(12) and (15). 
The s-wave alone produces a constant angular distribu- 
tion. The p-wave alone produces a distribution of the 
form 1 + a cos 2 9, i.e. symmetric around cos 9 = 0. For 
p-waves with the projection m = ±1 on the 10 Li mo- 
mentum, a is very small and negative, whereas a for 
to — is positive and of the order one. In other words 
the variation in the angular distribution is almost en- 
tirely due to the p-wave with projection 0. If both s 
and p-waves are present the distribution takes the form 
1 + a cos 2 9 + bcos9, i.e. becomes asymmetric due to 
mixing between these partial waves. 

Thus our pronounced asymmetry arises from the s and 
p mixing term, which dominates the angular variation, 
since the largest contribution from the s-wave alone is 
constant. A dominating p-wave would have produced 
a much more symmetric distribution. Substantial de- 
viations between measured and computed distributions 
would indicate a selective reaction mechanism empha- 
sizing specific partial waves. The shape of the angular 
distribution is independent of the beam energy as shown 
in the inset of Fig. 13, Only the absolute values of the 
cross section changes in accordance with Fig. & The 
shape in Fig. [l3| deviates slightly from that of pTJ due 
to the higher p-wave content. 
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FIG. 14. The invariant mass spectrum of 10 Li for the cases 
in Fig. ^. The rest mass is subtracted. The relatively small 
contributions from neutron and core participants are not in- 
cluded. The inset shows the distributions for different beam 
energies. 

The momentum distributions reveal properties of the 
initial three-body system and the reaction mechanism. 
The invariant mass spectrum of 10 Li after fragmentation 
of n Li on a light target carry in addition information 
about the properties of the two-body system. The com- 
puted spectrum is shown in Fig. [L4L where we only in- 
cluded contributions from the spectator neutron, i.e. the 



invariant mass E nc is constructed with the momentum of 
the spectator neutron even when the participant neutron 
is scattered by the target. In the center of mass system 
of 10 Li the invariant mass is, apart from the rest mass, 
simply the total kinetic energy of the two-body system. 

The spectrum must start from zero due to the phase 
space. The very low-lying peak is a signature of a domi- 
nating s-wave contribution. The height of the peak, not 
the position, reflects the energy of the underlying two- 
body virtual s-state. In contrast, a dominating p-wave 
contribution produces a peak at the energy of the cor- 
responding two-body resonance. In the present case the 
p-wave contribution is smeared out in a region around 
0.5 MeV. A narrow p-wave resonance would show up as 
a peak in this invariant mass spectrum [^2f . The mea- 
sured spectrum then indicates a low-lying s-state and a 
higher-lying p-resonance with a moderate width, see Fig. 
^c, where the curve computed with potential I is the same 
as shown in Fig. |IJ for a beam energy of 280 MeV/u. 
The contributions from stripping and diffraction as well 
as the energy dependence shown in the inset is consistent 
with the results in Fig. 0. The energy dependence of the 



peak heights and the widths are given in Table III 

As seen in Table |jll] all the cross sections first decrease 
with energy and then after a minimum slowly increase 
again. When the contribution from the scattered neutron 
is included the widths follow qualitatively the same pat- 
tern, although much less pronounced, and in particular 
they essentially do not vary for energies above around 200 
MeV/u. When only the spectator neutron contributes 
the widths are energy independent. This behavior dis- 
agrees with that of |20j], but is on the other hand consis- 
tent with the experimental results described in |T^,[l5| . 

In this comparison we implied that our FWHM is the 
same quantity as the width (or T-values) discussed in 
previous experimental papers. Although these quanti- 
ties are strongly correlated this is not strictly true, since 
the FWHM precisely is defined as the full width at half 
maximum of the computed distribution whereas the ex- 
perimental width is obtained as the width parameter in a 
fitted function. This easily shifts the emphasis from the 
small momenta in the calculations to the large momenta 
of the tail in the experimental analysis. This could easily 
produce uncertainties of several MeV in the FWHM. 

As we have demonstrated the distributions are more 
complicated than simple one-parameter functions. Com- 
parison of our FWHM with published experimental 
widths could then be rather misleading. The only safe 
procedure is to compare directly with the experimental 
data as we did throughout this paper. On the other hand 
then we immediately face the problem that the data too 
often include purely experimental effects related to beam 
profile, acceptance range, resolution and target thickness. 
It is therefore important to compare with the properly in- 
terpreted data or alternatively to fold the experimental 
effects with the calculated results. 
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V. DISCUSSION AND CONCLUSIONS 

Fragmentation reactions of weakly bound two-neutron 
halo nuclei provide detailed information about structure 
and reaction mechanism of corresponding three-body sys- 
tems. The three-body structure can be computed to the 
needed accuracy provided the interactions are specified 
whereas the reactions are much more difficult to treat 
propertly. Our aim is to understand and describe the 
principal features of all fragmentation observables qual- 
itatively and quantitatively to an accuracy comparable 
with that of the experimental data. As an important 
step we have chosen to investigate 11 Li considered as a 
three-body system in a physically simple and transparent 
model based on geometric properties and phenomenolog- 
ical interactions. 

One main assumption is often the sudden approxima- 
tion where the reaction is instantaneous. The reaction 
time must then be short compared to the intrinsic time 
scale of the relative motion of both the three halo parti- 
cles and the nucleons in the target and the core. Thus 
the intrinsic motion must be frozen during the collision 
or equivalently the beam energy must be large compared 
to the three-body binding energy and the Fermi energy of 
the target and the core. As expected this approximation 
has successfully passed the tests at high energies. 

In the present model we do not directly use the sudden 
approximation. We fully include the interaction between 
the target and each of the halo particles by use of a phe- 
nomenological optical model. The description of the in- 
teractions is then only limited by the validity range of the 
parameters employed in the optical model. On the other 
hand then only elastic scattering is described in details 
while all other processes are included as absorption from 
the elastic channel. Fortunately this is precisely the level 
of information required in discussions of fragmentation 
reactions, because the inelastic channels overwhelmingly 
produce different reaction products or particles scattered 
outside the detection range in the forward direction. 

Furthermore we only include the interaction between 
one halo particle (participant) and the target while ne- 
glecting the interactions between the other two halo par- 
ticles (spectators) and the target. In addition we also 
neglect the interactions between the spectators and the 
participant-target system in the final state. The halo par- 
ticles must then interact independently with the target 
as three spatially correlated but non-interacting clusters 
of nucleons. The motion of the spectators must remain 
undisturbed by the participant-target interaction. More 
precisely the two criteria for the validity of the model 
are that (i) the sum of the participant and target radii 
is less than the halo radius and (ii) the intrinsic veloc- 
ity of the partipant within the halo is much smaller than 
the velocity of the projectile or perhaps better that the 
characteristic time for the intrinsic halo motion is much 
larger than the collision time. 

These main approximations are justified for weakly 



bound and spatially extended halos colliding with a tar- 
get with an energy per nucleon larger than the intrin- 
sic kinetic energy of the participant within the halo, i.e. 
roughly the usual Fermi energy for nuclei although in 
principle unrelated to the Fermi energies of the nucle- 
ons within target and core. With these approximations 
one halo particle can interact with the target without 
disturbing the motion of the other two. The total cross 
section is obtained by adding the contributions from the 
three possible participants. The sudden approximation 
is reached when the optical model is reduced to a black 
sphere when the elastic scattering also is neglected. 

The model is described for point particles and the nec- 
essary generalization to finite radii involves the concept of 
shadowing. We eliminate the geometric configurations in 
the three-body wave function where the spectators move 
in the shadow of the participant. The need for this cor- 
rection arises from the simplifying choice of treating the 
participant-target interaction properly while leaving the 
spectators untouched in their initial state. If the spec- 
tators are close to the participant they would be either 
absorped or similarly scattered. However, these events 
contribute in the model computations with the probabil- 
ity given in the initial wave function. Consequently we 
must omit those unwanted configurations. Another im- 
proved treatment with better final state wave functions 
could directly take these effects into account. 

To use the model we must specify the interactions and 
the shadowing parameters. Within the halo projectile 
we have the neutron-neutron and neutron-core interac- 
tions supplemented by the three-body force. They are 
parametrized as gaussians to reproduce the 11 Li bind- 
ing energy, give a neutron- 9 Li ^ / 2 -resonance at 0.5 MeV 
with a width of 0.4 MeV and finally to produce about 
40% of ^-configurations in the n Li wave function. These 
requirements are necessary to reproduce various experi- 
mental data. The only freedom left for the halo inter- 
actions is then the spin-splitting, arising from the two 
different couplings to the core spin of 3/2, of both the 
S1/2 and pi/2-states in the neutron-core system. The 
related spin-splitting parameters influence neither the 
above data nor the fragmentation data. Good agreement 
with the data then indicates approximately correct sta- 
tistically averaged positions of the s and p-states in 10 Li, 
i.e. 0.71 MeV and 0.76 MeV. 

The two-body interactions between halo particles and 
target are described by use of the phenomenological opti- 
cal model with parameters adjusted to reproduce the cor- 
responding elastic scattering and absorption cross section 
data. The two shadowing parameters related to the sizes 
of halo particles and target are determined to reproduce 
both the absolute two-neutron removal cross section and 
the momentum distributions after fragmentation. 

We compute all momentum distributions related to 
fragmentation of n Li on carbon. For the same reactions 
we also compute the invariant mass of n Li and the neu- 
tron angular distribution, which recently was measured 
for 6 He fragmentation. These observables are in general 
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consistent with the available measurements. When the 
neutron and core participants are scattered they receive 
momentum tranfer perpendicular to the beam and the 
transverse momentum distributions are therefore broader 
than the corresponding longitudinal distributions. This 
is then a direct effect of the diffraction mechanism. 

The distribution for neutrons is relatively narrow due 
to the final state interaction, which affects the core less. 
The distribution for the center of mass of the core-plus- 
neutron spectators is the broadest reflecting the exten- 
sion of the initial wave function. Increasing the shadow- 
ing parameters decrease the widths of the distributions. 
The invariant mass reveals information about the low- 
lying continuum structure of the neutron-plus-core sys- 
tem. The large and very low-lying peak is the signature of 
a low-lying virtual s-state while the shoulder indicates a 
low- lying and fairly broad p-resonance. The angular cor- 
relation of the emitted neutrons in the neutron-plus-core 
center of mass system is highly asymmetric revealing that 
s and p relative neutron-core states roughly are equally 
populated in n Li. 

The experimental data are available for several ener- 
gies and targets, but systematic high precision data given 
as function of energy for one target does not exist at the 
moment. We can compute absolute values of a number of 
differential cross sections. However, in this paper we con- 
fined ourselves to the energy dependence of three-body 
observables for fragmentation on a carbon target with 
one excursion to an aluminum target. The distributions 
are essentially independent of target, but the absolute 
differential cross sections increase with target size. The 
scaling with target size seems to be somewhat larger than 
the square of the target radius. 

The energy dependence for the given carbon target is 
computed for a n Li beam of 50 to 900 MeV/u. The 
widths of the distributions are essentially constant above 
around 200 MeV/u and slightly increasing towards lower 
energies. The absolute values follow the participant- 
target cross sections. For neutrons this means a smooth 
increase with energy above around 150 MeV/u and a 
strong increasee towards lower energies. Furthermore 
diffraction contributes much less than absorption at en- 
ergies above 200 MeV/u whereas the inverse is true for 
energies below 50 MeV/u. This has the consequence that 
the widths of the transverse distributions are broader at 
low energies due to the domination of diffraction. 

These predictions presupposes that the model is valid 
at the low energies and the energy dependence of the pa- 
rameters are correctly included. The criteria for validity 
indicated relative neutron-target energies above around 
20 to 30 MeV and perhaps even lower for very pro- 
nounced halo systems. The optical model parameters 
for the neutron-carbon potential are adjusted to scatter- 
ing data down to these energies, but the shadowing pa- 
rameters are assumed to be constant. These parameters 
have a strong influence on the absolute values of the cross 
sections and a significant, but much less pronounced, in- 
fluence on the shapes of the distributions. The predicted 



widths could perhaps be systematically changed by small 
amounts due to such possible energy dependence. 

In conclusion, we computed systematically essentially 
all observables for the 11 Li three-body fragmentation on 
a carbon target. The same consistent model is used 
throughout. Most of the computations are in agreement 
with available measurements. This strongly indicates 
that the reaction mechanism essentially is correctly de- 
scribed in the model. The predictions are therefore use- 
ful as the unit for comparison with future experimental 
data. Detection of discrepancies would then be signifi- 
cant and therefore also very suggestive of necessary im- 
provements like, for eaxmple, different dependencies of 
some of the parameters. The same consistent model for 
all observables is crucial at the present level of accuracy 
and understanding. In this connection it is worth keep- 
ing in mind that treating 11 Li as a three-body system 
is an approxomation and the intrinsic structure must be 
unavoidable at some point in the interpretation. 



APPENDIX A: SPIN 1/2 SCATTERING ON SPIN 
ZERO TARGETS 

In Eq. (|^) we give the transition amplitude for the 
scattering process between the participant particle i and 
the target. Assuming that particle i has spin Sj = 1/2 
and the target has spin zero, we can write the transition 
amplitude as j33|. 

T^l = - — (Xs^\g(0) + h(8)(n-<T)\xs A ) , (Al) 

where 9 is the angle between p oi and p' Qi , er are the Pauli 
spin matrices and the functions g(9) and h(9) are given 
by Eqs. ([l8]) and (|l9|), respectively. The vector n is 
defined as 



n 



Pai x Poi 
\Poi x Poi\ 



(A2) 



(n ■ er) = (n • er)^ 



, (A3) 



which implies that (n ■ er) 2 = 1. 

We can then rewrite the key part of Eq. ( |Io| ) as 



E 



^2 



Mo* 



E (Xs^\g*{0) + h*{9){n-cr)\ Xs ^) 



*(Xs^>\g(e) + h{9)(n ■ ct)\ Xs ^>) 

x E K%»* M i 

Msj k T,j k 



fJM 
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(2tt) 



E^H^)l a + IM^)l a + 2W)^(fl)] 



SiE" 



(i^)i 2 + i^)i 2 ) e \ M > 



where we used that 

E (Xs,sJ(« ■ 0-)|X Si EJ' 



JM 1 2 



(A4) 



E;E" 



V M/ M s * E M/ M s E „ = 



(A5) 



as seen from Eq. ( |A3|) and the fact that the matrix 



is diagonal with identical diagonal elements j22|]. Inser- 
tion of Eq. @ into Eq. © and use of Eqs^l]) Jl|) 
and ( |l6| ) then immediately lead to Eq. ( |T^ ) and (13). 
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